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1. Introduction 



The c < 1 string theories, which can also be viewed as solvable models of 2D quantum 
gravity, have been solved perturbatively through their description as large matrix models 
(see the reviews Originally the matrix models were considered just as 'engins' that 

generate planar graphs, but very soon it became clear that the matrix models are in 
principle able to give a qualitative picture of the non-perturbative effectsS in these low- 
dimensional string theories 0. 

Recent works suggested to think of the matrix models as dual open string 

theories in presence of D-branes and led to an interpretation of the non-perturbative phe- 
nomena in terms of open string world sheets with appropriate boundary conditions. This 
interpretation was prepared, from the string theory side, by the remarkable works on Li- 
ouville CFT with boundary, which has been solved using conformal bootstrap methods 
2^26]. All comparisons between the matrix and CFT approaches showed agreement and 



led further to a new proposal for the matrix description of string theories with world sheet 
supersymmetry, type OA and OB theories p7| - p9| . 

In particular, the non-perturbative corrections (NPC) to the string partition function 
are believed to be produced by D-instantons [Q. On the string theory side, they were 
shown to be described in terms of a boundary CFT with Dirichlet boundary conditions on 
the Liouville field, known also as ZZ branes . Namely, the leading corrections are given 



by the exponents of the disk partition functions with ZZ boundary conditions []31|JT0| , pT | . 

For c < 1 minimal string theories, the leading non-perturbative effects were given 
a nice geometric interpretation in terms of a complex curve describing the closed string 



theory background [ p!2|p^ . Each point on this Riemann surface is associated with a FZZT 



brane, which implies Neumann boundary conditions for the Liouville field [^. The ZZ 
branes are associated with the double points of the Riemann surface, which can be also 
thought of as vanishing A-cycles. The disk partition functions with ZZ boundary conditions 
are given by line integrals along the corresponding S-cycles. 

The generalization of these results to c = 1 is not obvious because of the singular 
character of the limit c — 1 in the space of minimal string theories. In particular, the 
double points of the complex curve degenerate just to two singularities where all ZZ branes 
are situated. Therefore, in order to study the non-perturbative effects in c = 1 string 
theories, it is helpful to perturb the theory with a time-dependent tachyon potential, to 
the effect that the degeneracy of the singular points at c = 1 is lifted. 

The integrable tachyon perturbations are those with equidistant spectrum and it is 
natural to consider them in the context of a string theory with finite temperature. The 



^ The first systematic study of the non-perturbative effects in the simplest one-matrix model 
has been done by F. David El. 
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world sheet description of such a theory is given by the euchdean c = 1 string theory, 
defined by the world sheet action 



5c=i = — [ d'^ a \(dx? + (d(l))^ +2R(j) + ^(f>e^'^ + ghosts] , (1.1) 
An J I J 

in which the matter field x, or the local euclidean time coordinate of the string, is com- 
pactified: 

X + 27rR = x. (1.2) 

We are interested in the possible deformations of this theory by on-mass-shell tachyon 
operators^! 

Tp- /dVe^P^e(2-|pl)'^ (1.3) 



with the allowed by the compactification (1.2) values of the momentum 

Pk = k/R, /c = ±l,±2,... . (1.4) 
A general such deformation is achieved by adding to the action (1.1) a term 

dS = J2tk%/R. (1.5) 

We will be mainly interested in sine-Liouville deformation, which represents a perturbation 
by the two lowest vertex operators n = ±1. It contains, besides the Liouville term 
sine-Liouville interaction 

SS^^ = A Jd^a cos{x/R) e^^-Ti^'^, (1.6) 

where we denoted A = ti = t-i. 

The string theory deformed by the term (1.5) has been studied using the 'holographic' 
description provided by the Matrix Quantum Mechanics (MQM). The first exact result, 
the expression for the partition function for sine-Liouville deformation, was suggested by 



G. Moore in [^. Later this result was proved and generalised in [p3| - p5| using the fact 
that the perturbations (1.5) behave as To da integrable fiows. The classical background 
for given deformation is described by a complex curve whose exact form depends on the 
coupling constants in (1.5). This curve is not the one describing the FZZT brane, but is 



^ In the compact ified theory, in addition to this discrete spectrum of tachyon modes, there 
is a discrete set of the winding modes, or Kosterlitz-Thouless vortices. In the following we will 
consider only perturbations by tachyons; the results for perturbations by winding modes can be 
obtained by T-duality. 



2 



obtained from the latter by a projection that identifies an infinite number of sheets. The 
existence of two distinct complex curves is a peculiarity of the c = 1 string theory and is 
related with the logarithmic singularity of the resolvent at infinity. In the c < 1 string 
theories there is only one such curve. 

The leading non-perturbative effects in presence of sine-Liouville deformation were 
studied in [pil| , pT| and later, in |T^,|n[, where the string theory instantons were associated 
with the double points of of the complex curve. In this paper we study the quantum 
corrections to the leading NPC to the partition function, given by the pre-exponential 
factors. We calculate these factors using the matrix model description of the c — 1 string 
theory. 

We perform the calculation by two different methods. The first method uses the fact 
that the free energy of the perturbed theory as a function of the cosmological constant 
fx and sine-Liouville coupling A satisfies Toda partial differential equation. The second 
method is based on the formulation of the matrix model as a system of free fermions and 
consists in direct evaluation of the quasiclassical fermion wave functions. This method 
gives the answer for a general time-dependent background in terms of the complex curve 
associated with it. 

We are using the chiral formalism in which the role of canonical coordinate and mo- 
mentum are played by the left and right chiral combinations ~ x ±p. The advantage of 
this formalism is the exact bosonization of the fermion operators. Using the bosonization 
we interpret our results in terms of (target-space) CFT correlation functions. We show 
that the subleading corrections are related to the two-point correlation function of the 
bosonic field. Unlike the bosonization formulae for other matrix models for non-critical 
strings [^JT^,^, here we have two chiral fields associated with the left and right moving 
tachyons. The resulting bosonic field theory can be viewed also as a boundary CFT, in 
which the boundary condition relates the left and right fields, see e.g. [38,3^ for similar 
statement for the normal matrix model. Since our bosonic field describes chiral excitations 
above the Fermi sea, its correlation functions do not have direct interpretation in terms of 
the usual FZZT branes. Nevertheless, the one-point functions, evaluated at the singular 
points of the complex curve, reproduce the disk partition function on ZZ branes, as is the 
case in the c = matrix model WBM 



2. Time-dependent backgrounds in Matrix Quantum Mechanics 

2.1. The singlet sector of MQM as a system of free fermions 

The matrix quantum mechanics can be viewed as reduction of a two-dimensional U{N) 
gauge theory to one dimension, and involves one gauge field A^j and one scalar field X^j, 
both hermitian N x N matrices. The modern interpretation of the matrix path integral is 
as an effective open string theory on DO branes The matrix variable X describes 
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the open-string tachyon field in presence of the DO branes, and only effect of the gauge 
field is that it projects onto the singlet sector. The theory is described by the action 



5 = y"cftTr (PVAX-i(p2-X2)+^Ar), (2.1) 

where VaX = c^^X — i[A,X] is the covariant time derivative. The cosmological constant 
/i is introduced as a chemical potential for the size N of the matrices, which is treated as 
a dynamical variable. 

The closed c = 1 string theory appears as a theory of collective excitations of the 
matrix variables. The momentum modes arise as collective excitations of the matrix X^j, 
while the winding modes are collective excitations of A^j. The vertex operators (1.3) can 
be represented in the matrix model byi [^1|] 

fe-P*TrXlPl if»>0 
p |e-P*TrXH if p < 0, ^ ' 

where and X are the chiral combinations of the matrix coordinate and momentum 



X±P 

71" 



X± = (2.3) 



The theory simplifies significantly if formulated directly in terms of the chiral variables 
X_,_ |3^]. In the new variables the action (2.1) becomes 

S = j dtTr (X^ VaX_ + X^X_ + ^N) (2.4) 

so that the new Hamiltonian is linear in the canonical coordinates and momenta. 

In the singlet sector, characterized by absence of winding modes, the matrix model 
is described by a system of free fermions whose phase space is the spectral plane (x^, x_) 
of the two commuting matrices X^ and X_ . The fermions are governed by a one-body 
Hamiltonian 

Hq = —^{x_^_x_ + x_x_^_), (2.5) 



When comparing the results in CFT and matrix model descriptions, we encounter the old 



problem of operator mixing |4C]. This problem occurs because the correlation functions are 
integrated over the world sheet and the integrals have contributions from the coinciding points 
and a special prescription is needed to distinguish two tachyons with momenta pi and p2 close to 
each other from a single tachyon with momentum p\ + p2- This ambiguity possibly leads to an 
analytic redefinition of the couplings when passing to the matrix model description. 
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where 

£_] = —i. (2-6) 

The one-particle wave functions in "x^" and representations are related by Fourier 

transformation S: 

^_{xj = [^^J(xJ ^ rfx^e^^+--V,K). (2.7) 

The one-body Hamiltonian has continuous spectrum and is diagonalized, in the x_^_ 
and x_ representations, by the functions 

= x^'^-^ EeW. (2.8) 

The wave functions (2.8) have branch points at = and therefore are defined unam- 
biguously only for x^ > 0. It is useful to consider them either as multivalued meromorphic 
functions of the complex variables or as analytic functions of logx^. 

We will restrict ourselves to the case of a theory with one Fermi sea on the right of 
the top of the potential, in which the wave functions are supported, up to exponentially 
small terms, by the positive axis. Therefore we will define i/)'^ and i/'f by (2.8) along 
the positive axis, and their values along the negative axis will be obtained by analytic 
continuation from the upper and lower half plane, respectively.^ The wave functions (2.8) 
are orthonormal with respect to the scalar product 



{f,g)= / dx^fix^)gixj, (2.9) 





where the integration with respect to the chiral phase space coordinates is performed only 
along the positive real axis, x_^ > 0. The restriction of the phase space to > can 
be intuitively understood as follows. The wave function describes the system at given 
moment of time and its form for other moments is obtained by applying the evolution 
operator e**^. On the other hand, since the evolution operator shifts arg(x_,_) by =Fit, the 
analytic continuation of the wave functions (2.8) off the positive axis can be viewed as the 
result of evolution in imaginary time direction. 

The Fourier transformation S defined by (2.7) gives the reflection part of the scattering 
operator, relating the incoming leftmovers and outgoing rightmovers. It acts diagonally 
on the wave functions (2.8): 

S^E ^ e''^o{E)^E^ ^ ^2.10) 

with 

^^ME) = 1 g-f (E-i/2) YUE + 1/2). (2.11) 

The reflection phase 4>q{E), is the same as the one calculated using the standard x- 
representation of the upside-down oscillator in [^ |2[|. It has an exponentially small imag- 
inary part 

luiM-li) = ilog(l + e-2-^), (2.12) 
which determines the flow of particles under the potential barrier. 



In a theory with two Fermi seas one should introduce a second set wave functions, which are 
obtained from tp'^ by reflection Xj^ —Xj^. For details see the appendix of [35|. 
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2.2. Partition function and density of states 

The string theory compactified at time interval /3 = 2nR is described by the grand canonical 
ensemble of fermions at finite temperature 1/(3 and chemical potential /x. If is the 

partition function of the ensemble of fermions, then the free energy J-' — log Z is given by 

J^{fi)= [ dE p{E) log (l + e-^^^+^A , (2.13) 

J — oo 

where p{E) is the density of states related to the phase (j){E) of the fermion scattering by 

WA 1^ 

Zn Zn dE 

The non-trivial part of the density adds to a constant cut-off dependent term, where the 
cut-off A is introduced as the volume of the {x^ , x_ ) phase space. We remind the derivation 
of this relation in Appendix A, where we also recall the derivation of the relation between 
the free energy and the scattering phase 

2sm§^- J^ilj) = (P{-lj). (2.15) 

As we know from , this last identity allows to evaluate all tachyon correlation functions 
in the compactified theory, once their counterparts in the theory with R = oo are known. 

The identities (2.14) and (2.15) hold for any Fermi system at finite temperature and, 
in particular, after an arbitrary time-dependent perturbation that preserves the singlet 
sector. Thus all the information about the system is contained in the phase of fermionic 
scattering (t){E) as a function of the coupling constants {tk} associated with the tachyon 
operators (2.2). 

Let us recall the explicit expressions for the stationary background, where all t^ = 0. 
Then the phase (f}o{E) = (j){E)^tk=o} is given by (2.11) and the free energy itself is calculated 
by inverting the finite-difference operator in (2.15). The latter is diagonalized by the 
Fourier transformation 

where we reintroduced the cut-off A. Then (2.15) yieldJi 

•^^^^^^-°> = -4y. T sinhf sinh^ - ^^'^'^ 



^ Up to cut-off-dependent terms this integral gives the logarithm of the Barnes |43| Gamma 
function: J^(/i, 0) = \ogT2{^ -ip,\l,R). 
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The real part of this integral is the well-known integral representation for the partition 
function of the compact ified c = 1 string theory . Its genus expansion (which is also the 
expansion in l/^u) reads 



^pert(^){t,=0} = -f ^'log^ - ^log^ +i? J]^2-2h^;,(i?), (2.18) 

h=2 



where the genus h term Ch{R) is a known polynomial in 1/R. 

In the Fermi system, the tachyons are represented by collective excitations of the 
Fermi liquide P^J^]. A general time-dependent background, corresponding to the tachyon 
deformation (1.5) of the world-sheet CFT, can be introduced by changing the state of the 
fermi system in the infinite past and in the infinite future, in the spirit of the Polchinski's 
derivation of the tree-level tachyon S'-matrix. This amounts to modifying the asymptotics 
at infinity of the fermion eigenfunctions (2.8) in and x_ representations. The spectrum 
of the tachyons should be of the form (1.4) in order to have periodicity m. (3 — 2tiR in the 
imaginary time direction. As we know, such deformations of MQM are described by Toda 



hierarchy ||45| , |33| , |35[] . The commuting Toda fiows can be formulated either as an hierarchy 
of PDE with respect to the 'times' t±A;, or in terms of a pair of Lax operators satisfying a 
string equation Below we remind the necessary facts about the two approaches. 



2.3. Time- dependent backgrounds via Toda hierarchy 

The partition function of the perturbed theory is a r-function of Toda hierarchy and 
satisfies a hierarchy of PDE with respect to the couplings t±k. The first of them is the 
Toda equation, which is sufficient to describe sine-Liouville perturbation ti = t-i = A. 
Written for the free energy JF = log Z, this equation has the form 

i X-^dxXdxJ^ilJ, A) + exp [-4 sin^ (^|^) ^(^, A)] = 1 (2.19) 

and should be solved with initial condition provided by (2.18). 

Equation (2.19) defines the fiow between the critical points A = and /i = of 
the world-sheet CFT. When A is large, so that sine-Liouville term sets the scale, the 
cosmological term can be considered as a perturbation. In this case it is convenient to 
introduce the following variablesi 

y = fiC, ^={X^1^)-^. (2.20) 

The variable y is a dimensionless parameter and ^ ~ Qstr- Therefore, the genus expansion 
of the free energy is an expansion in ^ with y-dependent coefficients. For each coefficient 



This definition is adjusted to the case R < 1; the case R > 1 can be treated similarly. 
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in the expansion the PDF (2.19) yields an ordinary differential equation |3^]. Below we 
will only need the genus zero contribution jFg to the free energy. In the variables (2.20) it 
has the form 

^o(^,A) = ^ [§y^\og^ + fo{y)] , (2.21) 

the differential equation satisfied by /o can be integrated to an algebric equation for the 
function 

X(y)=c»J/o(y), (2.22) 

namely 

y = e-Ti^(y^ -e-^^^y\ (2.23) 
2.4- Time- dependent backgrounds via fermion wave functions 

As we mentioned, in the fermionic picture one can introduce sources for incoming and 
outgoing tachyons by changing the asymptotics of the one-particle wave functions. The 
perturbed one-fermion wave functions are obtained from the bare wave functions (2.8) by 
multiplying with a coordinate-dependent phase factor 

^Kx^) =e^*^±("±'^^i/'|(a;±). (2.24) 

The phases can be written as a sum of three terms 

^±{x^-E) = y±(xj + \(t^{E) + v±{x^-E), (2.25) 

where V± vanishes for = 0, v± vanishes for large x^ — oo, and (^(i?) is a constant. The 
deformations (1.5) of the world sheet CFT correspond, up to a possible analytic redefinition 
of the couplings, to the choice 

l^±(a:J = $^t±fcX^/«. (2.26) 

fc>i 

The constant mode </>(-£') and the pieces v±{x^]E) depend implicitly on t±n. They are 
determined from the requirement that the left and right fermions are related by the Fourier 
transformation (2.7): 

^^E^^B^ ^-l^fE=^B (2.27) 

where we absorbed the reflection phase factors e*"^*^^-* in the definition of the deformed 
wave functions. 

For large negative energy E, the compatibility of the saddle-point equations for the 
two integrals in (2.27) gives 

x^x_ = i 5^ x^J^ - E + j^Y.'^kiE) x;^/^ 

k>l k>l 

2.28 

= i 5^ x^_/« - E + i 5^ ^_.(i?) xl'^/^. 
fc>i k>i 
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This means that the two equations (2.28) define two functions = X^{x-^) which are 
inverse to each other: 

X^{X_{x^)) = x^. (2.29) 

The functions X^{x^), taken at the Fermi level Ep = — /x, can be considered as the 
gluing functions for an analytic curve Ai in C'^. This curve yields all the information 
for the given background. Its real section determines the shape of the Fermi sea in the 
quasiclassical limit n ^ oo. The functions X_^{x^) are in general multi- valued, but there 
is always a global parameter, the "proper time" variable r, such that the solution of (2.28) 
is given by x^ = x^{t). If all couplings t±k with k > kmax vanish, then the solution is of 
the form 

where the quantity x is related to the constant mode of the phase (2.25) via 

X = -R^<P. (2.31) 

The coefficients a±k and x can be found by substituting (2.30) into (2.28) and comparing 
coefficients in front of e r'^ . For example, in the case of sine-Liouville deformation, i.e. 
when only the first coupling constants are nonzero, this procedure gives 

/xe*^ - ^ (1 - A^e^^ = 1, a±i = A eT>^. (2.32) 

Taking into account that d^J-'o = x, it is easy to show that this equation is equivalent 
to the solution (2.23) of Toda equation. The solution beyond the tree level can be found 
using the full Toda integrable structure. 

The parametrization (2.30) can be thought of as a canonical transformation relating 
the phase space coordinates {x_,x^) and {t,E), where the proper time r appears as the 
variable canonically conjugated to the energy E: 

{x_,x^} = l ^ {t,E} = 1. (2.33) 

The origin of the last relation can be traced to the Lax representation of the operators 
x^ in the basis of the deformed wave functions, where the parameter uj = e"^ arises as the 
classical limit of the shift operator cD = e~*^^ p5|j^ . 
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3. NPC in sine-Liouville string theory from Toda equation 

In this section we restrict ourselves to the sine-Liouville deformation with coupling A, in 
which case the free energy can be determined from Toda equation (2.19). The initial 
condition is provided by the free energy of the non-deformed theory which is given by 
the integral (2.17). Before studying the theory with the sine-Liouville deformation, let us 
remind the exact expression for the NPC in the non-deformed c = 1 string theory following 
from the integral representation (2.17). The integral has a small imaginary part, which 
describes the flow of eigenvalues beyond the top of the inverse oscillator potential. It can 
be evaluated by extending the contour to the whole real axis and taking the residues at 
the two series of poles, = 2z7rn and = 2iTxRn: 

^op(^){t,=o} = t 4n(-l)"sin^ + ' ^ 4n(-l)- sin(7ri?n) ' ^^'^^ 



n 



This expression is of course compatible with (2.12) and (2.15). We see that there are two 
types of NPC, which have their origin in the two kinds of branes in the c = 1 string theory 
p!0[] . The flrst type is due to the D-instanton with Dirichlet boundary conditions on both 



matter and Liouville fields, whereas the second one is due to a DO-brane with Neumann 
boundary condition for the matter field. In both cases the pre-exponential factors do not 
depend on fx and therefore scale like g^^^.. 

The two types of NPC behave very differently with sine-Liouville coupling A. The 
NPC of the type e"^'^^"^^ do not fiow with A and are always given by the second term of 
(3.1). Indeed, given a solution of Toda equation (2.19), by adding a linear combination of 
exponents e~^'^^^^ with A-independent coefficients one obtaines another solution. 

In contrast, the NPC given at A = by the first term of (3.1), evolve non-trivially 



with A. The leading exponential contributions of this series were analyzed in [p!0| , pT[| for 
deformation by vortices. From these works we knowB that the leading non-perturbative 
effects have the form 

£^(^,A)~e-2«"(^)/^ (3.2) 

where 

9niy) = yOniy) + ^e-^ sin a = (3-3) 
and 9n{y) = dyQn is found as a solution of the algebraic equation 

sin^^= (i-1)"^ e^^(^^sin(i^^^). (3.4) 



^ The expressions obtained in [10,11| sliould be used after the T-duahty transformation ^ 



10 



The different solutions of (3.4) are labeled by the integer n in such way that 9^ — > 7rn 
when A ^ 0. 

Now we would like to find the subleading order, i.e. the factor in front of the expo- 
nential in (3.2). Following the procedure suggested in []T^, we take two solutions J-' and T 
of (2.19) that differ by the exponentially small quantity e = T — T . Then e must satisfy 
the linearized equation 

\\-^dx\dxe{ix. A) - 4 e-^^'^"^'^'^) sin^ (|^) £(^, A) = 0, (3.5) 

where in the exponent in the second term we approximated 

4i?2 sin^ (1^) A) ~ dlT^ = R\o^i + X{y). (3.6) 

This approximation is correct for our purpose since the subleading non-perturbative con- 
tribution is of order 0{g^t^ with respect to the leading one, whereas the perturbative genus 
expansion goes in powers of g1^^. Finally, changing the variables from (A, /u) to (^,|/) we 
write equation (3.5) as 

ocekydy + id^fe{i, y)=Ae-^ sin^ dy) e{C, y). (3.7) 

This equation is solved by a refinement of the Ansatz (3.2): 

£(^, y) ^ A{i, y) e-2^^(^)/^, A{i, y) = [Ca{y)]' , (3.8) 

where b may depend on y. We omit the index n since it does not appear explicitly in the 
equations. The details of the solution are presented in Appendix B. The result is that 
b = const and a{y) is given by (B.9). 

The constant 6, which gives the power of gstr in the pre-exponential factor, is not fixed 
by Toda equation. It looks like an integration constant. Therefore, the first idea is that it 
can be fixed from the initial condition (3.1) at A = 0, which implies that it should vanish. 
On the other hand, let us consider the limit of small A. In this limit y is large and 

9M-^n+{^-iy-' sm^y-"-^. (3.9) 

This leads to the asymptotics 

^(^, y) ~ i^/y)' y^ ~ A-^/^ (3.10) 

We observe that the subleading contribution does not have smooth small A limit. Therefore, 
the parameter b is not fixed by this approach. As we will see in the next section, it should 
be fixed to give the same power of ^'str as in c < 1 string theories, 

6=i, (3.11) 
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so that the full prefactor is given by the following function 



A{^,y) = C (e-^^ sin^^sin- 



) 



1/2 



2 e_ 

R 



[(i-l)cot(i^^)-cot^] 



(3.12) 



The overall coefficient C remains undetermined. Thus, we see that sine-Liouville deforma- 
tion changes drastically the behavior of the non-perturbative corrections in the subleading 
order: the power of Qstv ~ C changes discontinuously from to 1/2. This means that the 
limits A ^ oo and A ^ do not commute. As we will see later, the change of the be- 
havior of the non-perturbative corrections can be explained with the fact that the tachyon 
perturbations break the time translation symmetry of the theory. 

A nice consistency check of the result (3.12) is to reproduce the non-perturbative 
corrections to the free energy of the pure gravity in the limit where the couplings approach 
the c = critical point. We discuss this limit in Appendix C where we show that the result 
(3.12) does reproduce the correct critical behavior. 

4. NPC from the quasiclassical wave functions 

The approach based on Toda equation is conceptually straightforward but does not allow 
to fix the solution completely. Alternatively one can exploit the fermionic formulation 
to evaluate the non-perturbative corrections to the fermion reflection phase and then use 
the relation (2.15) to find those for the partition function. Moreover, this second approach 
allows to generalize the results obtained in the previous section to the case of a deformation 
with arbitrary number of non-vanishing couplings. Before considering the general case, we 
will illustrate the method for the stationary background, where the exact form of the 
non-perturbative corrections is known. 

4.1- The case of a stationary background 

Although the exact expression for the reflection phase is known, we follow here another 
approach which works only quasiclassically but can be generalized to other situations. For 
this we use the following matrix element of the scattering operator (2.10) 



where we used the orthonormality of the wave functions (2.8). Substituting the explicit 
expression for the wave functions (2.8) and introducing a cut-off A equal to the volume of 
the phase space, we write the diagonal matrix element as 



1 




(4.1) 



1 / / ax^ax_ ^ 




(4.2) 
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where po{E) Ri ■^\og{—K/E) is the density of states corresponding to the cut-off A. In 
the following discussion we can safely approximate po{E) by ^ log A. 

We would like to evaluate the integral (4.2) by the saddle point method. The two 
saddle point equations associated with the variables and x_ actually coincide and both 
give the equation for the classical fermion phase-space trajectory with energy E: 

x^x_ = -E. (4.3) 

The fact that the two equations coincide means that we do not have isolated saddle points 
but 'saddle contours' that go along the flat direction. If we change the variables to 

x±(T) = x/ie±^ (4.4) 

then the flat direction is along the 'proper time' r. The cut-off prescription < x_j_ < y/X 
then restricts the r-integration in (4.2) to the interval 

-ilogA<r<ilogA (4.5) 
and the integral in r gives the factor log — ~ 2iTpo{—e). The remaining integral in e is 



.A 

7o 



The saddle point is at e = —E and the gaussian fluctuations cancel \j2'Ke in the denomi- 
nator. This gives the genus-zero perturbative contribution for the reflection phase: 

0o(£;)f«£^log(-E)-£;. (4.7) 

The appearance of the non-perturbative corrections is related to the fact that the 
phase of the integrand in (4.2) is a multivalued function of e. As a consequence, there is 
an infinite number of saddles at 



log e„ = log(-£^) - 27rm (4.8) 

with n integer. The saddles with n > can be connected by constant phase contours to the 
dominant saddle and thus are potentially relevant. The integral along a contour that passes 
through the n-th subdominant saddle contributes a factor i^—\y^^'"f^E, ^%{e\o%{-e)-e) ^ 
where (—1)"^ comes from the pre-factor in (4.6). Taking into account the contribution 
of all saddles, one finds for the Fermi level E = —ji 

gi<^o(-M) g-i(MlogM-M) h + ^ g-27rnM 

\ neN 
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where the coefficients depend on the choice of the integration contour. This is in 
agreement with the expression that follows from the exact answer (2.12), 



-^) ^ -^log^ + ^ - i 5^ 2k (-1)" e-''^"^ 



(4.10) 



The coefficient ci = ^ in the contribution of the first subdominant saddle can be explained 
by the fact that the contour turns at 7r/2 after reaching the saddle |57]. The choice of the 
integration contours and the evaluation of the contributions of the subdominant saddles is 
a delicate problem and we will not discuss it here. 

Instead, we will try to understand the geometrical meaning of the saddles (4.8) in 
terms of the complex curve defined by the saddle point equation (4.3) at the Fermi level 
Ep = -II, 

x_^_x_ = fi. (4-11) 

As a real manifold the curve represents a hyperboloid, depicted in fig. la. It is characterized 
by two non-contractible cycles, the compact A-cycle, given by the section = x_, and 
the non-compact S-cycle given by the connected component of the section x_,_ = x_,_ with 
X, > 0, which coincides with the profile of the Fermi sea. The S-cycle connects the two 



infinite points oo_ = {x_ 
is given by the map 



00} with oo-i- = {x, = 00}. A global parametrization of Ado 



x^ 



(4.12) 



to the strip | Imr| < tt. 



a) 




contour 



X+X— — |-t 



physical 
sheet 



contour 
for X,. 



b) 




Rex 



Fig. 1: The Riemann surface of the non-deformed theory (one dimension is 
suppressed). The "saddle contours" are given by the i?-cycle and n copies of the A- 
cycle. The second picture represents the covering of the hyperboloid by the complex 
r-plane and non-minimal saddle contours for (r) in the universal cover. 
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However, because of the multi-valuedness of the integrand in (4.2), the solutions of the 
saddle point equation are actually described by the universal cover A4o of the curve (4.11), 
globally parameterized by the whole complex r-plane. For the minimal saddle, n — 0, the 
integration contours for dx^ and dx^ belong to the same sheet of the universal cover Mq 
and coincide with the cycle B on the hyperboloid, which is parameterized by the real axis 
in the r-plane. 

In contrast, the no n- minimal saddles, n > 1, are described by two different contours 
on the universal cover Aio: one of the contours describes x^{t) and the other one gives 
x_ (r). In the parametrization (4.12), they are represented as contours joining — oo — 27rm 
to oo and — oo to oo + 27rm, respectively, as the ones drown in fig. lb. These conditions 
ensure that Im(loga;_,_) = in the asymptotic regions of large The two contours 
have the same projection in Alo where they follow the cycle B and at some point wind 
n times around the cycle A, as it is shown on fig la. But since and x_ follow the 
different contours in the universal cover, the integrand in (4.2) along such a contour yields 
an additional factor {-l)^e~^'^'^^^. 

In this way, the non-minimal saddles are assiciated with 'double contours' on the 
hyperboloid (4.11), rather than with double points. The double points appear in the 
complex curve for the resolvent, y = w{x), which is also the c ^ 1 limit of the FZZT 
curves of minimal string theories. The curve for the resolvent is parametrized as 



and represents a Z2 orbifold of the universal cover of the hyperboloid (4.11), obtained by 
identifying the points r and — r. There are two points on the curve {±^/2i2, 0), which are 
also the positions of the branch points of y = w{x), which are images of infinitely many 
points on the r-plane r = mn. They appear as the (degenerate) limit of the double points 
of the minimal string curves. 

4.2. Quasiclassical wave functions in the general case 

Now we will apply the same method to calculate the NPC to the scattering phase in 
the case of a general tachyon deformation described in section 2.4. The first step of the 
calculation consists in the evaluation of the quasiclassical asymptotics of the fermion wave 
functions. 

First of all, we need the classical limit of the phases y:>±{x^,E) in (2.24), which is 
determined by the compatibility of two equations (2.28). The integration of these equations 
gives the following representation 





r sinhr 



(4.13) 



w 



77 



L 



X 



± 



^±{x^) 



X^{x'Jdx'^ +E\ogx^ +0±, 



(4.14) 



00 
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where the functions = X^{x-^) are the classical fermion trajectories defined in para- 
metric form by (2.30) and (p± are integration constants. As a result, in this approximation 
the wave functions take the form 

^|(a;±) ^S^(a;^)eT^'^±e^'i'I^^T^"±^'^"±, (4.15) 

where are factors coming from the subleading order in the Planck constant. These 
factors and the zero mode (p = 4'+ + 4'- can be fixed from the normalization condition 
similar to (4.1).i 

With the cut-off A introduced as in the previous subsection, the wave functions satisfy 
the normalization condition 

Va 



dx^^^ixj ^lixj = ^logA. (4.16) 







A second condition follows from the action of the scattering operator, that is the Fourier 
transformation. Since we absorbed the scattering phase into the wave functions, the Fourier 
transformation acts as the identity operator, which leads, together with (4.16), to 

dx^ / ^'s(2:_)e*^+^- ^'f:(x+) = p(E). (4.17) 

Jo Jo 

Substituting the asymptotic form (4.15) and keeping only the constant cut-off dependent 
piece of the density (2.14), one writes the two conditions as 

dx^\Bj^ = ^logA, (4.18) 
dx^ I c^a;_S_S+e*^^"+'"-^ = -^logA, (4.19) 



'0 JO 

where we introduced the effective action 



5(x+,x_) = x^x„ - / A:„(x;)dx; - / X^{x'_)dx'_. (4.20) 

The integral in the second condition can be evaluated by the saddle point method. 
The property (2.29), following from the compatibility of two equations (2.28), means that 
the leading contribution to the integral is again associated with a one-dimensional saddle 
contour 7f in the {x^,x_) plane, defined by the functions x^ = X_^{x^) and going from 
X = oo to X = oo. For E = —fi the saddle contour defines the profile of the Fermi sea. 



Note that the deformed wave functions are not orthogonal anymore. 
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To evaluate the gaussian fluctuations in the transversal direction to the saddle contour, 
it is convenient to change variables from {x^,x_) to (r, — e), where r is the same as in 
(2.30) and the variable e parameterizes the transversal direction. The equation of the 
saddle contour in the new variables is e = —E. For the quadratic form 6'^S{x^,x_) one 
finds, using (2.33), 

5^,Six^,xJ =(2d,x^d,x^ - ^ {d.x^Y-^ (d.x^Y) ^ 

" ('4 211 

d^x^ d-rX_ J 2 d-rX j^d-rX _ 2 

As a result, the condition (4.19) is written as an integral along the saddle contour 7f: 



g'J^p-^-"-+ / dr B_B^ [-drx^drx_y'^ = '^\ogK. (4.22) 
The two conditions (4.18) and (4.22) are satisfied if 

(/> = (/>pert(£^) = / x„ B^ = \ . (4.23) 

In particular, taking into account the relation (2.15), this result gives the expression of the 
//-derivative of the perturbative free energy J^^ert as an integral over the S-cycle on the 
complex curve (2.30) Knowing the wave functions (4.15), one can also reconstruct the 



icj) 



quasiclassical wave function in the x-representation. We refer to Appendix D for details. 



4.3. Instantons as double points of the complex curve 

In order to find the non-perturbative corrections for the free energy it is sufficient to 
calculate those for the scattering phase 4>{E). This can be done by taking into account 
also the non-minimal saddles of the integral (4.19). We saw in section 4.1 that if there is no 
tachyon potential, the non-minimal saddles are given by pairs of one-dimensional contours 
living on different sheets of the universal cover of the hyperboloid x^x_ = /i, but having 
the same image in C^. The collective coordinate along these saddles is the Minkowski 
time, since the background is time-independent. In presence of tachyon potential there is 
no time translational invariance. As a consequence, the minimal saddle contour will evolve 
with time and the non-minimal saddles will be given by isolated saddle points, or double 
points of the deformed complex curve, to be described below. The meaning of these double 
points is essentially the same as in the minimal string theories ||T^ , p^ . 

Again, the saddle-point equations x^ = X^(x_) and x_ = X_{x^) define a complex 
curve Ai in C^, which is a deformation of the universal cover of the hyperboloid x^x_ = 11. 
The Riemann surfaces of the functions and X are the projections of the curve M. to 
the x_ and complex planes. Due to the multi-valuedness of the functions X^{x_) and 
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X_{x_^), equation (2.29) has also particular solutions describing isolated double points, 
where the complex curve touches itself [1T3[| . 

The double points can be classified most easily by using the uniformization parameter 
r, defined in equation (2.30). They represent pairs of points r' 7^ r" in the complex r-plane 
such that 

x^{r') =x^{t") and x^{t') = x_{t"). (4.24) 

Each double point is by construction a solution of (2.29) and as such represents a saddle 
point for the action (4.20). 

Let us restrict ourselves to the case = t-k for all k, which is simpler to analyze. Then 
Ofc = tt-k in (2.30) and the set of double points is given by the pairs (r' = —iOn, r" = z6'^), 
where Q^, are determined by 

h 

sin {Q^ = ^ ttk sin 0^) , 9n ^ ttu when aU 0, (4.25) 

fc=i 

As t' and r" are pure imaginary and complex conjugated, x^{±i9n) = x_ {±i9n) are both 
real. For the sine-Liouville deformation the discrete set of parameters 9n is determined by 
the equation 

sin(^„) =a(y)sin(i^0„) (4.26) 

with a{y) from (2.32). This equation is the same as (3.4) and therefore the parameters 9n 
coincide with those defined in section 3. 

Let us calculate the contribution of the n-th saddle point to the normalization integral 
(4.19). The leading contribution is given by the value of the action (4.20), which depends 
on the integration contours for dx^ and dx_ along the complex curve. The two contours 
can be parameterized in the r-plane by the intervals 

r G (00, 0) U (0, -z^n) forx,(r), 

+ (4.27) 
r G (-00,0) U (0,z6'„) forx_(r). 

As a result, the action is given by S'pert + iSn, 

Spert = X dx^, Sn = i (p x_ dx^, (4.28) 

Jl¥ Jin 

where the contours 7f and 7^ are respectively the images of the intervals (— cxd,cxd) and 
(i^n, —'i'9n) in the r-plane. 

The first term Spert gives, according to eq. (4.23), the leading perturbative approxi- 
mation to the scattering phase. Its integration contour is the S-cycle connecting the two 
punctures at x^ = 00. On the other hand, the image of the integration contour 7^ for the 
second term Sn is a closed loop in C^, because it ends at the n-th double point. If the 
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double point is considered as a vanishing A-cycle of the complex curve, then the contour 
is a the dual compact S-cycle. 
To find the subleading pre-exponential factors, one should evaluate the fiuctuations 
around the saddle points. The variation of the action (4.20) gives 



S'^Sn{x , x_) = 5x,5x_ 



dx 



2 dx 



(5x^y 



dx 



2 dx 



iB. 



{5x_f. 



(4.29) 



This leads to the following prefactor in the calculation of the integral (4.19) 



dxj^ 

~9r 



dx _ 



i9„ 



1 - 



dx 



dx 



-1/2 



(4.30) 



where the first two factors come from the quasiclassical wave functions. 

Adding together contributions to the left hand side of (4.19) of the minimal saddle 
contour, given in (4.22), and all isolated saddle points, we geti 



oo I — 



Sn 



e^'^(-^) log A. 



(4.31) 



Taking the logarithm of both sides, one gets a quasiclassical expression for the zero mode 
that includes the non-perturbative corrections 



x_ dx^ — i 



7F 



E 



TvAr 



n> 



^0 v^logA 



(4.32) 



Finally, we can apply (2.15) to evaluate the non-perturbative terms for the free energy 
itself which gives 

^ - -^Pert + . if . y: e-'- . (4.33) 



2V21ogA;^„sin%|^ 
The classical action Sn associated with the n-th double point can be written as 



Sr. 



E 

X d-rX^dr = — 2 / O^de, 



(4.34) 



where we used the canonical transformation (2.33), and the prefactor to the exponent has 
the form 



An ~ 



sm 



2R 



sm 



2 e_ 

R 



iOr. 



dr 



-i9„ 



-1/2 



(4.35) 



^ Here we do not specify the contours of integration and therefore the combinatorial factors in 
front of the individual terms. We assume that the contribution of each saddle point enters with 
the factor 1/2. 
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The free energy is determined by (2.15) up to terms of the form e~'^'^^^. The analysis 
of section 3 based on Toda equation shows that these corrections do not depend on the 
tachyon potential and therefore are given by the last term on the r.h.s. of (3.1). However 
there is an argument, advanced in [|1^, that such terms do not appear whenever a tachyon 



potential is switched on. 

The result (4.35) gives the instanton corrections for an arbitrary tachyon deformation 
of the theory. Now let us restrict ourselves to the case of sine-Liouville deformation. 
It was already shown in [|l3l that in this case Sn reproduce the leading non-perturbative 
corrections (3.2) obtained by Toda equation. Therefore let us concentrate on the subleading 
contribution (4.35). Calculating the derivatives using (2.30), where only terms with k = 1 
are present, and taking into account the defining equation (4.26) for 9 = 9^ one finds 

r 1 1"'/^ 
A{^,y) = C le-R^ sin^ 9 sin^ {{^-l)cot{^ 9) -cot 9)\ , (4.36) 

where the overall coefficient C is found to be 

C= 'f^ . (4.37) 
4v^logA 

This result is in complete agreement with the result (3.12) obtained by integrating Toda 
equation. 

In particular, this confirms our claim that for any finite tachyon deformation the 

1 /2 

subleading corrections scale like g^l,. and thus they are non-analytic in the limit A ^ 0. 
This is due to the breakdown of the time translation invariance in any time-dependent 
background. As a consequence, in the first case the contributions to the saddle point 
approximation come from one-dimensional saddle contours, whereas in the second case 
they arise from isolated saddle points. Hence, the determinants of fiuctuations, giving the 
main non-trivial contribution to the subleading correction, are one- and two-dimensional, 
correspondingly. This explains the difference in the power of the string coupling. 



5. The instanton effects from bosonization 

In the matrix model, the tachyon modes are collective excitations of fermions propagating 
as left and right moving waves on the Fermi surface. After second quantization, these 
modes form a bosonic field. In {x,p) representation of MQM, the fermions are non- 
relativistic and as a consequence the bosonic field is self-interacting. In contrast, as in 
(x_^,x_) representation the one-particle fermionic Hamiltonian becomes first order, the 
fermions can be exactly bosonized. This bosonization is essentially the one that occurs in 
Toda hierarchy, modulo some subtleties related to the fact that we are using a non-compact 
realization of the latter. 

In this section we show that the non-perturbative corrections to the free energy can 
be expressed in terms of vertex operators for a boson field. The bosonic field in question 
describes not the whole spectrum of tachyons in the theory, but only the discrete subset 
that survives after the compactification. The bosonic field formalism works particularly 
well in the quasiclassical limit and gives a nice interpretation to the results obtained in the 
previous two sections. It might also help to interpret these results in terms of D-branes in 
c = 1 string theory with time-dependent background. 
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5.1. Second quantized fermions and density operator 

We start by introducing the second quantized fermion fields 



(5.1) 



where the operator amphtudes b(E) and b^{E) satisfy the canonical anti-commutation 
relations 

{bHE),b{E')} = d{E-E'), (5.2) 

with all other anticommutators equal to zero. We work with only one set of operator 
amplitudes, since the left and right fermion operators are related by Fourier transformation 
and the fermion wave functions have the property (2.27). The last property followed from 
the fact that we absorbed the reflection phase in the definition of the wave functions, so 
that the fermion reflection operator acts as the identity operator on the amplitudes in 
^'-representation. 

We are interested in operators in the matrix model of the form 

Of= Tr/(X^,X_), (5.3) 

where /(x^, x_) is a smooth function of its two variables. The operator Of is well defined 
as the matrices and X commute due to the gauge field. In the second quantized 
formalism, this operator translates into 



where 



J dx^dx_f{x^,x_)U{x_^_,x_), (5.4) 



U{x^,x_,t) = ^e'^+^-^'t_(2:_,t)^'^(x^,t) (5.5) 

is the fermion phase space density operator. The expectation value of the operator W, or 
the Wigner's function, is evaluated with respect to the thermal vacuum of the compactified 
theory which is defined by 

{^\bHE)biE')\^) = ^f-=^. (5.6) 

In particular, for the trace of the identity operator N = {fx\Tr l\fj,) we find, using (5.6) and 
(4.17), 

N = j dx^dx_W{x^,x_)\^i) = ^ • (5.7) 
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Comparing with (2.13), one reproduces the weU known relation between the number of 
fermions and the grand canonical free energy 



where we explicitly included the non-universal cut-off dependent term. Then (2.15) allows 
also to write the relation of to the reflection phase 

sin {^d^) N = (log A + a^(/>(-^)) . (5.9) 

5.2. Bosonization formula for the fermion wave functions and the density operator 

The fermionic operators can be expressed as exponents of a bosonic field with continuous 
spectrum of energies. It is however technically more advantageous to introduce another 
bosonic field with discrete spectrum associated with the possible energies in the Euclidean 
compactified theory. If we restrict the spectrum of tachyons to be discrete as in the 
compactified theory, the spectrum of the bosonic fields will be also restricted to the discrete 
set of purely imaginary momenta pn = in/ R, ti G Z. The sector of the theory spanned on 
these states gives a 'non-compact' realization of Toda integrable structure.0 The bosonic 
field is thus spanned on the operators {E, Oe, t±n, dt_^^} and has the form 

^JxJ = V±{xJ + ^4>-E\ogx^+D^{xJ, (5.10) 

where V± are the potentials (2.26) and the differential operators (p and are defined by 

^ = -idE. ^±(^j = E^^r/''^*±"- (5.11) 

n>l 

The deformed fermion wave function at level E = — /i, known in the mathematical 
literature as Baker- Akhieser function, can be written as the expectation value of normal 
ordered exponential of a bosonic field. A straightforward generalization of the bosonization 
formula for the compact Toda hierarchy gives BTI 



^-M(x^) = (27rx±)-^/2 . gT*-l±(x±) . |^\ ^ ^512) 



Usually in Toda hierarchy the wave functions are bi-orthogonal polynomials and the r- 
functions are labeled by a non- negative integer s £ Z-|-, the degree of the polynomial [49,50|. 



In our case the role of s is played by the discrete complex variable ^ + i{n + 2)"r' ^ ^ 

use the terms compact and non-compact Toda hierarchy in the analogy with the representations 

of the compact su{2) and the non-compact s/(2,M). 
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where the normal product sign : : means that all derivatives are moved to the right, and 
the expectation value of the differential operator O is defined as {O) = • O • Z. From 
the bosonic representation of the Becher-Akhieser function we get the following operator 
formula for the fermion bilinears 

= (^: e-'*-(^-^)^ e-**+(^+^, (5.13) 

where the hermitian conjugation changes the normal to anti-normal ordering in the first 
factor. Passing to the normal ordering and subtracting a divergent term in the exponent, 
we finally obtain 

^ 2h{x+x_)-^ •e-^^K'^'-) (5.14) 
where we introduced the full bosonic field 

^{x_^_,x_) = ^_^_{x_^_) + ^_{x_). (5.15) 

Since we performed a subtraction, the overall coefficient on the r.h.s. of (5.14) was fixed 
by hand. This is done by comparing the quasiclassical expressions obtained below through 
bosonization with those that follow from the quasiclassical wave functions. 

Using (5.14), one can easily write the bosonization formula for any observable of the 
form (5.4). In particular, for the trace of the identity operator (5.7), which gives the 
number of particles, one obtains 

5.3. Quasiclassical limit 

The quasiclassical asymptotics of the fermionic wave functions is determined, through the 
representation (5.12), by the quasiclassical expansion of the exponentials of the bosonic 
fields $_,_(a;_^). Here we will restrict ourselves to the first two orders: 

where 

^JxJ^i^JxJ) (5.18) 

is the vacuum expectation value of the field $_,_(x_,_) and (■ ■ ■)c denotes the connected 
correlator. In the leading order the phase of the fermion wave functions are thus given by 
the expectation values of the left and right bosonic fields: 

(fi^{x^;-li) + li\ogx^^^^{x^). (5.19) 
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In particular, the quasiclassical phase of the fermion scattering is equal to the expectation 
value of the zero mode cj) 

cl>={4>). (5.20) 

In order to find the subleading factor in the exponential, we need the two-point cor- 
relation functions. The connected two-point correlators are generating functions for the 
derivatives d^d^J-' and dndj^J-', which were calculated in the case of sine-Liouville defor- 



mation in from Toda hierarchy. The answer is actually true for any deformation and 
can be expressed in terms of the functions t(x^) obtained by inverting the parametric 
representation (2.30) [0.0 Note that the functions t(x_,_) are given by the derivative in 
/U of the phase (5.19) of the fermion wave functions 

t(x J = ±a^$ J. (5.21) 
The explicit expressions of the generating functions are 
(x J =-^dlT- log (x^ e^^(^±)) , 
D^{x^)b^{y^)T= ^dlr + \og{x^-y^)-\og (^e±^ - e±^^ , (5.22) 
D^{x^)D_{yJJ^= \ogil-e— — 

From the first two identities (5.22) one finds the connected correlation function for the left 
and right chiral fields: 

(: I.,(xj|.,(,j :) = log L- ^^-^p7(^) + 4^5^ (5.23) 

V 2(x±y±)2Hsmh — ^ ^ / 

Substituting this in (5.17), one obtains the quasiclassical form of the fermion wave func- 
tions: 

^-^{xj = , (5.24) 

Up to the last term in the exponent, the expression (5.24) reproduces the previously 
obtained result (4.15) together with (4.23). 

To explain the appearance of the last term, let us consider the zero modes 4>± . To all 
orders, they are given by 

eT^<^± = (eT*^/2) ^ ^(^^ 27^) ^ (5 25) 



In the case of a compact Toda hierarchy, when the fermion eigenfunctions are entire functions, 
these relations have been obtained in |52,^|. 
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so that for the total zero mode we reproduce the expression (2.15): 

^i4> ^ gi(<^++<^-) ^ ^{i^+Jr) (5 26) 

On the other hand, the quasiclassical calculation gives 

^Ti4>± eTi(0/2)-|(0<^)^ ^ (5_27) 

which is, of course, agrees with the expansion of (5.25). In this way the last term appears as 
a correction to the zero modes of the wave functions of chiral fermions. This term reflects 
the freedom in the choice of the relative normalization of the left and right fermions and 
is cancelled in the product 

Now let us turn to the fermion bilinear (5.14). In the leading order, its expectation 
value is given by the exponential of the expectation value of the full bosonic field, 

{^{x^,x_)) ^ ^^{x^) + ^_{x_). (5.28) 

To find the subleading contribution, one should again evaluate the two-point correlation 
function of ^{x_^_,x_). From (5.22) one finds 

. 1 T(,r:_)-r(:(/ ) . T{y _)-t(x_) 

{:^(x.,x )^(y.,y ):> =log , , , , , . 

\ ^ ^^-^-^ /e %inhl^f±^^sinh^^^4^ (5.29) 

+ iog(x;i/^ - 2/;i/«) + iog(x:i/« - 2/:^/^), 

which yields for the expectation value of the fermion bilinear 

(//|*t(a;_)#4xj|//) = ^(a;^a;_)-^e-^*^"+'"-)-^(^^^"+'"-)^^"+'"-)^)c 

g-i*(a;^,x_) ^ (5.30) 

^ 2^^;^;^ 2i? sinh ' 

The last two terms on the r.h.s. of (5.29) appear due to the normal ordering. Without 
normal ordering the correlation function of the gaussian field would be given only by the 
first term. Taking this into account, the formula (5.30) generalizes straightforwardly to 
the case of several fermion bilinears: 

(/.I n mxi)^4xiM = n ' det ] I . (5.31) 

/=! L\ 2n^drx1drx'^_ ^'^ V 2i2 siuh ^^^^^ ' 
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We expect that this formula becomes exact if one replaces the function u± = e^^^±' by 
the shift operator e*^^ of the Toda hierarchy.lll 

Note that the second factor in the denominator of (5.30) diverges on the surface 
of the Fermi sea where t(x_^) = r[x_). This divergence appears due to breakdown of 
the quasiclassical approximation of the bosonization formulae near the Fermi sea. This 
divergent factor can be canceled by application of a difference operator in ii. Indeed, in 
the given approximation such an operator acts only on the exponent. Therefore, applying 
the relation (5.21), one concludes that 

2Rs\n{^d^) Ui>l{x_)i>^{x^)\ii) = - \ I ~J =. (5.32) 

\ / 2nW—drX,drX_ 



As a result, the r.h.s. becomes the quasiclassical expression for the product of the wave 
functions \1/I'^(x_ )\l/^'^(x^). This allows to write a difference equation for the expectation 
value of any operator Of of the form (5.3) 

2Rsm{^d,){^\d\^^) = -j^ I ^7^^^/K,^Je^^+^--^''(^+'^-). (5.33) 

The non-perturbative effects to the free energy follow from this result if one specify 
O to be the identity operator (/ = 1). Then the relation (5.9) implies 

dx,dx_ — -j^^^=^^^^ = log A. (5.34) 



^^'^ J y^—drX^drX 

The evaluation of this integral by the saddle point method gives the non-perturbative 
corrections to the zero mode (j) and, through the relation (2.15), to the free energy itself, 
as it was done in section 4.3. 



6. Discussion 

In this paper we studied non-perturbative corrections to the partition function of the 
compactified c = 1 string theory deformed by a generic tachyon source. The flows between 
different backgrounds are described by a non-compact realization of Toda hierarchy. The 
space of all such backgrounds is parametrized by a set of Toda 'times' {t±fc}, which in 
the world sheet CFT have the meaning of coupling constants for the allowed marginal 
deformations. Each such time-dependent background is described by a complex curve, 
which gives the solution of Toda hierarchy in the dispersionless, or quasiclassical, limit. 
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In the case of the two-matrix model, where the integrable structure is that of KP hierarchy, 



such an exact formula was derived in |21]. 
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We studied the dual reahzation of the string theory as the singlet sector of the matrix 
quantum mechanics, described by free fermions in upside-down oscillator potential, using 
chiral canonical coordinates ~ x±p. In the fermionic system the complex curve appears 
as the complexified classical trajectory at the Fermi level. This curve has an infinite (for 
R irrational) number of double points, obtained as solutions of a transcendental equation. 
The leading non-perturbative corrections are associated with these double points, in the 
same way as in the minimal string theories. The exponents for the leading NPC are given 
by integral of a holomorphic differential along the compact cycles associated with these 
double points. 

Our main new result is the expression for the subleading correction given in (4.35) for 
a general tachyon potential and in (4.36) for the particular case of sine-Liouville theory. We 
derived these results by quasiclassical analysis of the fermionic wave functions as well as 
by solving the linearized Toda equation in the simplest case of sine-Liouville deformation. 
Then we showed that these results have a natural interpretation in terms of a bosonic 
collective field, whose oscillator modes are in correspondence with the Toda couplings tk- 
The fermion bilinears then can be represented as vertex operators of this bosonic field, 
eq. (5.14). The subleading corrections are proportional to the exponent of the two-point 
bosonic correlation function, eq. (5.30). 

It is natural to expect that the NPC have their origin in the D-branes of the string 
theory. In order to have such an interpretation, they should be expressed in terms of world 
sheets with local boundary conditions. Such an interpretation exists for the eigenvalue x of 
the random matrix variable, which is the boundary cosmological constant for the Liouville 



On the other hand, our results were derived in terms of the chiral phase-space coor- 
dinates, and X , which do not have direct interpretation in terms of local conformal 
boundary conditions in the world-sheet theory. Therefore, to find the meaning of the non- 
perturbative corrections as the effects due to D-branes, we should first translate the results 
obtained in the (x_^,x_) representation back to the (x,p) representation. 

Assuming that tk = t-k and that there is a finite number of non- vanishing couplings, 
the parametric form of the curve in the {x,p) space is obtained directly from (2.30): 



k=i 

The FZZT curve is the one for the resolvent y = w{x). The latter is related to the 
momentum p by 



field. 




k=l 



(6.1) 




p{x) 




(6.2) 



2i 
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This equation is easy to solve only in the case of stationary background, where the two 
sides of the cut of w{x) are the images of the straight lines Imr = ±7r, and whose solution 
is given by (4.13). In the case of a time-dependent background the two sides of the cut 
of w{x) are parametrized by the lines T±{t) = t ± Z7ri(t), t > 0, where 7ri(t) tt only 
asymptotically at t —>■ ±00. Therefore eq. (6.2) can be written more explicitly as 

^^^^ ^ w{t + iTViir)) - w{t - iTVijT)) 

where the function ni{t) is obtained as the solution of the transcendental equation 

sinh ffl — — )t) 

sinTTi + Vafc ^ sin ((1 - |)7ri) = 0, t eM. (6.4) 

^-^ smht ^ " 

k=l 

As a consequence, the parameter r does not uniformize the curve y = w{x). In general, 
w{t) will have infinitely many branch points. 

Considered at t = 0, eq. (6.4) coincides with the condition that drX = for purely 
imaginary r. Therefore the point r = ±z7ri(0) is a branch point for p{x). It is easy to 
see that all purely imaginary solutions of (6.4) correspond to branch points of p{x). The 
solutions of (6.4) can be classified by their asymptotics at infinity: 7r„(t) mr. 

The equation (6.4) for 7r^(r) is not compatible with the transcendental equation (4.25) 
for the parameters 9n of the double points, except for the stationary background. There- 
fore, for time-dependent backgrounds the double points never occur at the branch points 
of the resolvent or along its cuts. This fact can be geometrically understood by consider- 
ing the classical trajectory x{i9) in imaginary time direction. The branch points are the 
turning points of the classical trajectory, while the double points are associated with its 
intersection points. Therefore, turning on a time-dependent tachyon background resolve 
the infinite degeneracy of the double points of the complex curve for the c = 1 string 
theory. 

Now let us return to the interpretation of the non-perturbative corrections in terms of 
D-branes. For the leading non-perturbative effects, which are given by the closed contour 
integrals Sn (4.28), this was done in |T^,|n[. Using the relation (D.IO) of the light-cone 
coordinates to the spectral density and (6.2), one obtains 

Sr, = -\i dx {w{e'^x) - wie-'^'x)) . (6.5) 

The resolvent w{x) is equal to the derivative of the FZZT disk partition function with 
respect to the boundary cosmological constant x: 

w{x) = d^^^^^^ix). (6.6) 
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Taking into account the symmetry ti;(e*^x(r)) = w{e *^x(— r)) of the resolvent, this aUows 
to identify 

-S^ = ^^,,^{e-'^x{ien)) - $^^^^(e^"x(z^„)). (6.7) 



In [|T3[ it was also checked in the first order in the deformation coupling A that the 



leading correction is given by the ZZ disk partition function 

-S^ = ^,,{n,l) (6.8) 

so that the n-th double point of the complex curve is associated to the (n, 1) ZZ brane. 
Assuming that this equation remains to be true in the full deformed theory, from (6.7) and 
(6.8) we conclude that the ZZ partition function is equal to the difference of two FZZT 
partition functions calculated from two sides of the double point. This is a generalization 
to the case of a general tachyon deformation of a similar identity in the non-deformed 
theory (A = 0). iHSUl • 

Note that our analysis also reveals the distinguished role of the (1,1) ZZ brane. Up 
to now it was the only one which played a role in the interpretation of matrix models as 
theories of open strings. In our context it is associated with the first double point of the 
complex curve. One can show that this is the only double point which is situated on the 
physical sheet of the Riemann surface. It is natural to expect that this fact corresponds 
to the observation in CFT that the (n, 1) ZZ branes with n > 1 possess some pathological 
properties [^. However, if we are calculating the non-perturbative corrections to observ- 



ables that explore lower sheets of the Riemann surface, than the other double points may 
become relevant. 

The world-sheet interpretation of the subleading order of NPC seems to be a more 
complicated problem. It is expected to be related to the annulus amplitude on ZZ brane. 
However, the expression given for this amplitude in |T^] diverges in this case. One may 
attempt to relate this divergence to the factor (log A) ^ in (4.37). But in this way we can 
not explain the universal non-trivial functional dependence on the couplings. On the other 
hand, our results are naturally expressed in terms of the two-point amplitude of a gaussian 
field whose left and right components are associated with the incoming and outgoing 
tachyons. It is not however obvious how to reformulate them in terms of the annulus open 
string amplitude on FZZT branes, which is defined in terms of the eigenvalue x. 
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Appendix A. Relations between the density of states, the scattering phase and 
the free energy 

Knowing the zero mode (^{E) of the fermion wave functions, one can reconstruct the density 
of states by introducing a completely reflecting cut-off wall at distance A ^ // . The 

wall introduces a boundary condition at = a;_ = 

[^^'](VA) = ^'(VA). (A.l) 

Thus, one identifles the scattered state with the initial one at the wall. Applying this 
condition to the deformed wave functions (2.24), one obtains that it is satisfied for a 
discrete set of energies En (n G 2) defined by 

(PiEn)-En\ogA + V{A) + 2nn = 0, y(A) = ^^(tfc + t_fc)A'=/2«. (A.2) 

k 

From (A.2) one can find the density of the energy levels in the confined system 

2n 2n aE 

Now the free energy J^{fj,,R) = log Z{fj,, R) of the fermionic system compactified at 
distance (3 = 2%R can be calculated using (2.13) with the density (A. 3). Dropping out the 
A-dependent non-universal contribution and integrating by parts, one obtains 

H^) = -^l dct>{E) log (l + e-^(^+^^) = -R dE^^^^j^. (A.4) 

We close the contour of integration in the upper half plane and take the integral as a sum 
of residues. This gives for the free energy 



^(/x)=^ <Pi^r/R-^i) (A.5) 

from which (2.15) follows. 



r=n+i>0 



Appendix B. Solution of (3.7) 

In this appendix we solve the equation (3.7) using the ansatz (3.8). Substituting (3.8) into 
(3.7), keeping only the leading terms in the ^ — > limit (which come from the terms with 
all or all except one derivatives acting on the exponent in (3.8)), and taking into account 
that g{y) satisfies the equation 

(1 - ydy)g{y) = sin Hdyg{y)] , (B.l) 
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one finds the following first order differential equation on A{^,y): 



a [2 ((1 - ydy)g) {ydy + ^d^) log A-{1- ydyfg] 
= -i sin ^ log^ - ^ cos ^. 

Taking into account the form of A{^,y), one obtains two equations for a(y) and 



(B.2) 



h' 



2ae1fi{l-ydy)g + ^sm^ 



a [2b ((1 - ydy)g) (1 + ydy log a) - (1 - ydyfg] = -^sm^ dyloga - ^ cos ^. 

(B.3) 

The only solution of the first equation is 



b = const. 



Then the solution of the second equation is given by the following integral 



(B.4) 



logo 



(.) = / 



dy 



ae^{l- ydyfg - ^ cos ^1 - 2a (1 - ydy)g 



(B.5) 



2a eT^ y{l - ydy)g + sin ^ 
Using (B.l), the integral can be rewritten as 

{±{2+^dyX)-l) 



log a{y) = -^log ae^y{l - ydy)g + ^sin^ 



+ 



/ ^^1 F 

J ^ cos % 



(B.6) 

Then we change the integration variable to X using the equation (2.23). The second term 
in (B.6) becomes 



dX 



2R-1 
bR 



cosh(^X)+ (-^-l)e^^-e-^ 



-X 



J cos ^ - 2Vasinh {^^X) 



(B.7) 



Changing the variables to 6 = g' by means of (3.4), after simple algebra one arrives at the 
following integral 



-/ 



sin 6* sin I sin (i^^) 



log 



C sin I 



sin6'sin(i^6') 



2R , 

H log 

2i? - 1 ^ 



sin^ 



L(sin|)^-Msin(l-«^))« 



4-1 



(B.8) 
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where C is an integration constant. Combining this result with the first term in (B.6) 
rewritten in terms of 6, one obtains 



log a(y) = ^ - 2 log 



sin I 



- ^log [{^ - 1) cot {^e) - cote] -\\og [C-^sin^sin(i^^)] . 
As a result, the pre-exponential factor reads 



(B.9) 



2<_ 



sin^ I 



{sin^ ^sin^ e){{i- 1) cot 6) - cotO)}-'/' . 

(B.IO) 



Appendix C. The c = critical limit for the subleading contribution 

It is well known that sine-Liouville theory exhibits a critical behavior when sine- 
Liouville coupling becomes sufficiently large. Physically, at the critical point the fluctua- 
tions of the matter field get frozen in minima of sine potential. As a result, the system 
approaches the c = CFT describing pure two-dimensional gravity. 

In it was shown that the leading non-perturbative correction to the partition 

function corresponding to n = 1 follows the same pattern, whereas other corrections with 
n > 1 disappear in the c = critical limit. Thus, we expect that Ai[y) must exhibit a 
singularity a,s y ^ yc with 

yc = -{2R-1)R~^{1- R)^ (C.l) 

and, furthermore, the behavior of Ai near this singularity should reproduce the non- 
perturbative effects of the c — Q theory 

Ai-glll^.-iy-yc)-'''. (C.2) 

Indeed, the singularity ai y = yc corresponds to a critical point of (2.23), near which 
the relation between y and X degenerates: 

^ ^ - + 0{{X- X^f) . (C.3) 

yc 2R 

At the critical point ^ and the first two terms in the expansion of 9i around the 
singularity are 

e^{y) = V3{X^ - XY'^ - "^^^^^o" f ^ (^c - Xf'^ + O ((X. - Xf'^) . (C.4) 

Due to this result, the leading term in (3.12) is 

rM^ ((I - R)X)~2£^ 

Ai{y) ^ C \^ , ' ' {X, - X)-''/\ (C.5) 

Taking into account (C.3), one finds the law (C.2) corresponding to the c = theory. 
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Appendix D. Quasiclassical wave function in the ^-representation 

In this appendix we study the quasiclassical wave function of the deformed ferniionic system 
in the ^-representation where x is the usual fermion coordinate related to eigenvalues of 
MQM. The simplest way to get this function is to use the result (4.15) together with (4.23) 
for the same wave function in the light-cone representation. Then it is enough to apply a 
unitary operator relating the two representations. 

This unitary operator can be represented as an integral operator with kernel defined 
as a solution of the following equations 

—^{x^idx)\x){x_^\=x^\x){x^\, (D.l) 
v2 

It is easy to check that it is given by 

W{x-J = ^e^*('"-^^"±+<). (D.3) 
Then one obtains (the same result will be for the "x " representation) 



'^^{x)= / dx^\x){x^\'^1{x^) 
Jo 



(D.4) 



The integral can be calculated by the saddle point method. The saddle point equation 
reads 

a;^ +X_(a; ) x^ -^-(a^+) ,^ 

— vi — = — T2 — = ^ 

where X (a;_^) is the function defined by (2.30). The fiuctuations of the phase around this 
pomt are — -|- — f — • Altogether these results give 

*-(a;) = ^-1— , exp ii [ p{x') dx') , (D.6) 
^^2n\^rX\ \ J J 

where the low limit of integration in the phase can be chosen to cancel the zero mode 0+. 

Finally, let us show how the relation, which is evident in the classical limit, between the 
density q{x) of eigenvalues and the functions x^ = Xj^x^) follows from the quasiclassical 
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wave function (D.6). The density in the compactified theory can be written in terms of 
the fermion wave functions as follows 



,w=2 (D.7) 

— oo 

where the coefficient 2 comes from the fact that there are left and right moving fermions. 
Substituting the quasiclassical limit (D.6) for the wave functions and replacing the thermal 
factor by 6'(— /x — E), one finds 

d^g{x) = (D.8) 

7r\OrX\ 

From here, using the identities 

dx. 



dE 



one gets the well known result for the quasiclassical spectral density 

q{x) = :;^{x^-x_), x= ^{x^+x_) . (D.IO) 
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